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Abstract 



We study a unitary version of the one-dimensional Anderson model, given by a five 
diagonal deterministic unitary operator multiplicatively perturbed by a random phase ma- 
trix. We fully characterize positivity and vanishing of the Lyapunov exponent for this 
model throughout the spectrum and for arbitrary distributions of the random phases. This 
includes Bernoulli distributions, where in certain cases a finite number of critical spectral 
values, with vanishing Lyapunov exponent, exists. We establish similar results for a unitary 
version of the random dimer model. 

1 Introduction 

Unitary operators arise naturally in quantum mechanics as the time evolution of the Hamil- 
tonian in solving the time dependent Schrodinger equation. In particular, for Hamiltonians 
with periodically time-dependent potentials, the spectral properties of the monodromy op- 
erator (the unitary operator giving the evolution over one time period) are the central 
object of mathematical investigations (e.g. [T], [9], [H]). A model using random unitary 
operators was used in [3] to study the single particle behavior of an electron in a small 
one-dimensional metal ring in the presence of a large uniform electric field generated by 
a linearly ramped magnetic flux. The unitary operator used in describing this system is 
characterized by a five diagonal band structure. Motivated by this model, the spectral 
analysis of a class of random unitary operators with similar band structure was undertaken 
in [U, [IS], [E] and [l3]. 

This class of operators on /^(Z) can be written (up to a unitary equivalence) in the form 



where S is a unitary operator with a five-diagonal matrix representation. S depends on a 
parameter t G (0, 1) which controls the size of its off-diagonal elements and takes the role of 
a disorder parameter for Uuj, see Section[2]below. is a diagonal matrix of random phases, 
diag{e~*^fc }. For our application, {6*^ : A; G Z} is a sequence of i.i.d. random variables on the 
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one-dimensional torus T = M/27rZ with a non-trivial probability distribution fi, i.e. supp/i 
has two or more elements. U^j can be considered as a "unitary Anderson-type model", 
where S plays the role of the free Laplacian and where the perturbation is introduced via 
multiplication rather than addition to ensure that the resulting operator is still unitary. 
Indeed, due to the band structure of the operator, the generalized eigenvectors can be 
studied using complex 2x2 transfer matrices. This formalism allows to introduce the 
Lyapunov exponent 7(A), where e*'^ is the spectral parameter, see equation ([8]) below. Due 
to the fact that the transfer matrices have determinants of unit modulus, the Lyapunov 
exponent is almost surely non-negative. This also allowed to prove a unitary version of 
the Ishii-Pastur Theorem, i.e. Theorem 12.11 below . and to deduce the absence of absolutely 
continuous spectrum in case of uniform distribution /x [1]. In [H] it was shown that the 
spectrum of f/^ is almost surely pure point with exponentially decaying eigenfunctions if the 
distribution ^ has a non-trivial absolutely continuous component and non-empty interior. 
This holds for arbitrary value of the disorder parameter t. Thus, up to this point, the 
results for one-dimensional unitary Anderson models are in close analogy to those known 
for the self-adjoint Anderson model. 

Further results on unitary Anderson models were found in [16j, where the density of 
states is studied, and in [18], which develops a fractional moment approach to prove local- 
ization for multi-dimensional analogs of unitary Anderson models. 

The main goal of this paper is to further investigate the model ([1]) and to fully determine 
the set {A G T : 7(A) > 0} for arbitrary disorder t and arbitrary distribution fi, including 
singular and, in particular, Bernoulli distributions. 

For most choices of fi we find that 7(A) > for all "quasi-energies" A € T, see Theo- 
rem [3]2j However, there is one exceptional situation: If is a Bernoulli measure supported 
on two diametrically opposed points, i.e. supp/i = {a,b}, \a — b\ = n, then there exist 
two critical quasi-energies A = —a and A = —b at which the Lyapunov exponent vanishes, 
while it is positive for all other values of A (Theorem 13. ip . In fact, we show in the proof 
of Theorem 13.11 that at the anomalies A = —a and X = —b the transfer matrices Tn{uj,X) 
satisfy the asymptotics 

-E((ln||r„,(a;,A)||)2) ^00, (2) 
n 

i.e., roughly, ||r„(c^, A)|| ~ e^^")'^'. 

As there are no more than two critical values of A, the unitary version of the Ishii- 
Pastur Theorem shows that the unitary Anderson model ([1]) almost surely has no absolutely 
continuous spectrum, irrespective of the underlying probability measure and disorder. We 
expect that methods such as those used in [8j can be adjusted to also show that C^j almost 
surely has pure point spectrum with exponentially decaying eigenfunctions, i.e. is spectrally 
localized, but we haven't carried out the details of this. 

It is interesting that the structure of Lyapunov exponents for the unitary Anderson 
model is richer than for the self-adjoint one-dimensional Anderson model. For the latter 
it has been long known that the Lyapunov exponent is positive at all energies for all non- 
trivial single site distributions of the random potential, e.g. [5]. 

In a more general class of self-adjoint Anderson-type models it has been shown that 
the existence of critical energies with vanishing Lyapunov exponents can lead to the co- 
existence of spectral localization and suitable forms of dynamical delocalization, e.g. |15j . 
The simplest self-adjoint model which shows this phenomenon is the so-called dimer model 
|12tlll|. in which the random phases appear in the form of identical neighboring pairs. The 
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typical anomalies encountered in the dimer model are stronger than those in Theorem 13.11 
below in the sense that transfer matrices, rather than satisfying (jj]), are uniformly bounded 
in n. In [15] it is shown that this leads to super-diffusive transport, while, to our knowledge, 
the dynamical effects caused by an anomaly as in ([2]) have not been studied (they should 
be much weaker, if detectable at all). 

In Section [6] we study a unitary version of the dimer model, where we can show that the 
Lyapunov exponent is positive away from at most finitely many critical values. However, 
for the dimer model with Bernoulli distributed phases, i.e. supp/x = {a,b}, and such that 
\a — b\ is in the spectrum of S, there are two critical values where transfer matrices are of 
the type studied in [15], in particular they are bounded in n. 

While we do not carry out a complete study of spectral and dynamical localization 
properties of unitary Anderson models, in Section [7| we state one more result which has 
been relevant in this context in the self-adjoint case, namely continuity of the Lyapunov 
exponent in A away from the critical quasi-energies. This is proven by a rather direct 
adaptation of the proof in the self-adjoint case, e.g. [5]. 

Let us finally mention that unitary operators with the same band structure as S and 
Uuj above also arise in the form of so-called CMV-matrices in the study of orthogonal poly- 
nomials on the unit circle, e.g. [19\ [6]. In this setting, a definition of Lyapunov exponents 
different from ours arises naturally (applicable for example for the case of i.i.d. random 
Verblunsky coefficients), see Section 10.5 in [19]. However, it can be shown that both 
definitions lead to the same value [17] (if either one of the two Lyapunov exponents exists). 



2 The Model 



Analogous to the self-adjoint case, we look at a random unitary operator as a random per- 
turbation of a deterministic ("free") unitary operator. The model and the results presented 
in this section can be found in [1], [13], [16] and [18|, see there for details and additional 
motivation and results. Motivated by [3], we choose the free unitary operator S on P{Z,) 
with band structure 
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where the position of the origin in Z is fixed by (e2fc_2, 5*62^) = — t^, with {k G Z) 
denoting the canonical basis vectors in /■^(Z). The real parameters t and r are linked by 
r'^+t'^ = 1 to ensure unitarity. Due to unitary equivalence it suffices to consider < t,r < 1. 
Thus S is determined by t. We shall sometimes write S{t) to emphasize this dependence. 
Excluding trivial special cases, we assume < t < 1. The spectrum of S{t) is purely 
absolutely continuous and is given by the arc 

a{S{t)) = S(t) = {e'^ : i9 G [-arccos(l - 2*^), arccos(l - 2^^)]}, 

which is symmetric about the real axis and grows from the single point {1} for t = to 
the entire unit circle for t = 1. 
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The random perturbation is then introduced via multipHcation by a diagonal matrix 

Z)^ = diag{e-*^^}, (4) 

with {0^ : /c G Z} a sequence of i.i.d. random variables on the torus T = ]R/27rZ. More 
precisely, we introduce the probability space (O, J^, P), where is identified with T^, T is 
the (T-algebra generated by cylinders of Borel sets, and P = (2)fcGZ '^here is a probability 
measure on T. The random variables Qk on (17, T ^ P) are defined by 

6lfc : ^ T, e% = 0Jk, Z. (5) 

This ensures that the resulting operator 

= D^S (6) 

is unitary and ergodic with respect to the 2-shift in Q [4]. U^^ also inherits the band 
structure of the original operator S and has the almost sure spectrum [16] 



E = exp(i supp ^)S(t) = {e*"S(t) j a G supp 
Here supp /i denotes the support of the probability measure defined as 
supp// := {a I /i(a — e, o + e) > for all e > 0}. 
Solutions of the eigenvalue equation 

Uuii) = e'^ip, = CkCk, 

with Cfc G C, A E C, are characterized by the relations 
for all /c G Z, where the transfer matrices T : — > GL{2, C) are defined by 

and the phases by 

et{x) = et + x. 

Note that det T{9^^{X),0^i^^^{X)) = e*(^2fc-e2fc+i) has modulus one and is independent of A. 
We have for any n G N 

=n^2(n-i)(A),^2Vim(^))---^(^o(A),^r(A)) -t„ka) 

=nr2„(A),r2„+i(A))-i...r(r2(A),ri(A))-i ^T_„(a;,A) ('-i 

We also set Tq{uj, A) = I. 
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As shown in [3], for any A G C, the Lyapunov exponent 

7^(A)= hm ^ln||r„(a;,A)|| 

n^itoo |n| 

almost surely exists, has the same value for A; — > oo and fc — > — oo, and takes the determin- 
istic value 

A version of the Ishii-Pastur theorem suited to the present model was proven in 
Theorem 2.1. Let U^j he defined by ^ and ^ and 7(A) by Then 

SacC{e^AG5i;7(A) = 0} . 

All norms on GL{n,C) being equivalent, we choose to work with the row-sum norm 
for convenience. Thus, in what follows the norm is the maximum row sum, i.e. for A = 
(^«i)i!j=i' ll^ll ~ niaxi<j<„ ^"^-^ \aij\. P(C^) denotes the projective space of C^, we write 
V G P(C^) for the direction of w G \ {0}. The action of a 2 x 2 matrix A on P(C^) is 
defined by Av = Av. 



3 The Main Results 

For a particular choice of the underlying distribution of the random phases the unitary 
Anderson model, unlike the self-adjoint one, exhibits two critical values of the spectral 
parameter, where the Lyapunov exponent vanishes. 

Theorem 3.1. // suppfj, = {a, b} and \a — b\ = tt, then 

(i) 7(-a) = li-b) = 0. 

(ii) 7(A) > 0, for aUXeT\ {-a, -b}. 

This is the only exceptional case. For all other choices of the probability measure // the 
Lyapunov exponent never vanishes. 

Theorem 3.2. // {a, b} C supp fj, such that \a — b\ ^ {0, vr}, then for every X €T we have 
7(A) > 0. In particular, if supp jjL contains at least three elements, then 7(A) > for all 
A G T. 

Theorem 13.21 is a generalization of the corresponding results on the Lyapunov exponent 
previously proven in [3] and [13] . 

Theorem I3.ir i) will be proven in Section SI while the proofs of Theorem I3.1f iil and 
Theorem 13.21 are provided in Section [3 

The fact that for all non-trivial probability measures /x, the set of critical quasi-energies 
contains at most two points combined with Theorem 12.11 gives the following immediate 
corollary, concerning the almost sure absolutely continuous spectrum of U^. 

Corollary 3.1. For any non-trivial distribution fi of the i.i.d. random phases, we have 
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4 Critical Quasi-energies 

In this section we prove Theorem 13. l( i). i.e. that a Bernoulh measure /i with diametrically 
opposed masses a, b indeed gives rise to two critical quasi-energies at A = —a, A = —b. 
Denote /x(o) = p £ {0, 1) and = q = 1 — p. 

For A = —a, the i.i.d. random matrices 7'(^2A;(-^)) ^2A:+i('^)) ^^^^ °^^y following 
values with non-zero probabilities, 

'r(0, 0) = -/, with probability p^ 

(1 2r/t \ 

„ , with probability q^ 
2r/t (3r2 + l)/ty 

T{iT, ~ ( Q 1 ) ' with probability pq ^ ^ 

T(0, vr) = ( j , with probability pq. 

2r/t —1 J 

The latter matrices take much simpler forms when represented with respect to the basis 
i{ f \ ^ /+ ) ' ( I- ^1 \ /+ ) } °f 'C^- Hence, we define the matrices A{9, rj) as 



(r + l)/t;' V(r-1)A 



It follows that 

^(0,0) = -/, A(7r,7r) 



P 

1/p 



with p := (r -h l)^/*^ > i. 

A straightforward calculation shows that 

, E('"l|A,.(^.A)ll)^ 



where An(w,A) = U'^^^A{9'^^{X) , 

^2fc+i(A))- order to simplify the notation, we will 
suppress the to dependence of various quantities for the remainder of the section. 

v^n,2(A)^ 



Let no := and Un{X) 



:= A„(A)mo- 



Lemma 4.1. If Xn = — ^ ^'^^ for n > 0, i/ien 

In p 

7(-a) = Inp lim -E{\xn\). (12) 

n^oo n 



Proof. From (jlOp it follows that with probability one there are x, y G M with xy = 1 and 
either A„(— a) = ( g ^ 1 oi' A„(— a) = ( g 1 • both cases it follows readily that 

||An(-a)|| = ||A„(-a)no||oo, 
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where || • ||oo denotes the max-norm on C^. This imphes that 

7(-a)= hm iE(ln||A„(-a)no||oo). (13) 

n— +00 fi 



Furthermore we see from the specific form of A„(— a) that A„(— a)no 
Therefore 

1 

In ||A„(-a)no||oo = h-imax(|n„,i(-a)| 



Un,i{-a) 
l/un,i{-a))' 



Un,i{-a)\ 
= 1 hi \unA{-a)\\. 

The required result then follows from (jl3p and the definition oi Xn- □ 

The following lemma is devoted to the necessary analysis of the random sequence x„,. 

Lemma 4.2. {xn)n>o is an integer-valued Markov chain with xo = and transition prob- 
abilities 

F{Xn+l = Xn) = P^, = Xn + I) = , 

P(x„+i = -Xn) = pq, F(Xn+l = - iXn + 1)) = PQ- (14) 

Proof. Clearly, xq = 0. Let A{n + 1, — a) := ^(^2(n+i) ~ ^2{n+i)+i ~ o,))- Iii the case 
A{n + 1, —a) = —I, we have |?x„4.i^i(— a)| = |u„^i(— a)|, i.e. Xn+i = Xn- If A{n + 1, —a) = 
p 



Q / ' ^^'^^ \un+i,i{-a)\ = p\un,i{-a)\ and Xn+i = x„ + 1. Similarly, A{n + 1, -a) 



1 

1 



implies that |?i„,+i^i(— a)| = \un,2i—o-)\ = l/|n„,^i(— a)| and thus Xn+i = —x^ 
-l/p\ . I / M _ 1, / M _ 1 



Finally, A{n + l,-a) = ( ^ \^^] gives |u„+ii(-o)| = -|'u„,2(-a)| - , ^ 

\-p U y p p\un,i[-a)\ 

i.e. Xn+i = —Xn — 1- Thus Xn+i is determined by x„ and A{n + 1, —a). The transition 

probabilities follow from Q and p^ . □ 



(15) 



Lemma 4.3. n — > oo, 

E(x^) ^ q 
n 2p 

In particular, we have that for all n, 
n\xn\) < Cn^^^. 

Proof. Let a = q — p, we denote by ¥,{x\y) the conditional expectation of x given y. It 
follows that 

E{Xn\Xn-l) = 0?Xn-l + qOl. 

Since IE(a;„) = E(E(x„|x„_i)), we get that 
E(x„) = a^E(x„_i) + qa 
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and, iterating, 



1 - a^" 
E{xn) = a^"E(xo) + qa- 



Similarly, since 

E{xl\xn-i) = xl_i + 2qxn^i + q, 
we have that 

E{xl) = + 2qE{xn^i) + q. 

Another induction gives that, for all n > 3 



l-a^'^^, , ,^ 2qa 



E{xi) = E{x'o)+2q- E{xo)+nq{l + —^: 

1 — 1 — 

+ 2q'^a[l - —^{2 + a^ 



1 1 _ "-2) 

o[2 + a — 5 — 

I - 1 - a2 • 



Since \a\ < 1 and Efxo) = IE(xq) = and qi\ H — ^-l = (7/2p, we get (fTSll . Using that 

1 — 

E(|x„|) < (E(x2))i/2^ in turn, proves the second assertion and finishes the proof. □ 

Lemma [4. 3 1 is. in fact, a consequence of general extensions of the Central Limit Theorem 
used in the study of dynamical systems, e.g. Section A. 4 of We include the previous 
elementary proof for the convenience of the reader. 

The main result of this section now follows immediately. 

Proof of Theorem \3.1^ ). The fact that 7(— a) = follows directly from (fT2]) and 
Lemma [43l The proof of 7(— 6) = is identical. □ 



5 Positivity of the Lyapunov Exponent 

In this section we show that, except for the two critical energies discussed above, the 
Lyapunov exponent ^ is positive. This constitutes the contents of part of Theorem 13. II 
and of Theorem 13.21 For each A G T, the random variables 9q and 9i induce a measure 
on GL(2,C) through T{e^ + A,6l^ + A). Denote the smallest closed subgroup of GL(2,C) 
generated by the support of this measure by Ga,^. Thus Ga,^ is generated by the matrices 
T{9, rj), defined in ([7]), where 9 and rj vary in A + supp /i. 

Fiirstenberg's Theorem [2] states that if Ga,/^ is non-compact and strongly irreducible, 
then 

n— >oo n 

The proof that Ga,/^ is non-compact for all values of A G T was given in [T3] and holds 
for any non-trivial probability distribution fi. For completeness we repeat the proof here. 

Lemma 5.1. Ga,^ is non-compact. 
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Proof. Let 9 and t] be on the torus, ^ t], and let x := e z := e Let G{d, rf) be the 
closed group generated by T{9,9), T{ri,ri), T{9,r]), and T(t],6). Define 

D:=Ti9,9)Ti9,rj)-'=(^^^^__^^ 1^gG{9,v), (16) 
E := T{rj, 9)-'T{e, 0) = ( J |f " "^^^ ) G G(0, rj), (17) 

^:=^^=f^f- 1, ^^^';.|^ ) GG(0,r?), (18) 

\ ^(xz — 1) xz — ^Ixz — 1| / 

J .- ED = ( f~ ^^"i^ ~ '^^^ "^^M G G(0,r/). (19) 

\ j(l — xz) xz J 

Note that detL = det J = 1 and that J'^ = L*. Thus we get the self- adjoint element 
K := J~^L of G{9,rj). In fact, K is positive definite and detiC = 1. More calculation 
shows that 



2r^ ,2 

tri^ = l + ^xz-ir + 
= 2 + —\xz-l\^. 



r 



2 



xz 7t\xZ — 1 



t4 

As 9 r] and therefore xz 7^ 1 we conclude that trK>2. Positivity of K implies that 
it has an eigenvalue strictly bigger than 1. Thus, containing all powers of K, the group 
G{9, rf) is non-compact. Li particular, with 9 = \ + a and r/ = A + 6, we see that G^,^ is 
non compact. □ 

It remains to prove strong irreducibility under the assumptions of Theorem I3.1l![iil) as 
well as under those of Theorem 13.21 Under the already established non-compactness of 
Ga,^, strong irreducibility of G\^^ is equivalent to 

#{5^ : 5 G Ga,^} > 3 for ah G P(C2), (20) 

see [2]. We first use this fact to prove that for supp^ = {a, h} and \a — h\ = vr, —a and —h 
are the only critical quasi-energies. 

Proof of TheoremlEMM}- Let 9 = a + X, r] = b + X, with A G T\{-a, -b}. In the 
terminology introduced above, the condition that |a — i| = tt can be written 
{— 1, 1}. Since x = — the operator L G Gx^^ defined by (fTHj) takes the form 



L 



-1 
-2r 




As detL = 1 and [trL| > 2, L is hyperbolic, hence iterations of L map any direction in 
P(C2) to infinitely many directions, except when v coincides with the direction of one of 

its eigenvectors, given by t>+ = r + 1 , w_ = r — 1 . 

Next we prove that even for the eigenvectors of L, we have that ^{gv^ : g G Ga,/^} > 3. 
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Under the current conditions, the transfer matrices take the form 



Tiv,e) = 

Therefore, we have that 
/ 




T{e,e)v+ 



T{rj,0)v+ 



r + 1 



(r — x) 



r + 1, r(r + 1)^ 

-[rx + X) + 



t3 
r+1 



t3 



r + 1 



[rx 




x) + 



r + 1 



(r + x) 



\ 



r + 1 . _. r 
-—^[rx - x) + - 



A simple calculation shows that T(6,6)v+ = T{9,ri)v-^ only if x = r. Similarly, 



T{9, r])v = T{ri, 9)v only if x € {—1, 1}, while assuming that T{9, 9)vj^ = T{ri, is equiv- 
alent to (r + l)t'^ = 0. All these cases are excluded by the assumptions of Theorem [SJjju]) . 
Therefore, we conclude that T{9,9)v+, T{9,r])v+, and T{ri,9)v-^- are all different. In a 
similar way one treats V-. We thus have proven that 

#{gv:geGx,^,}>3 for all G P(C2). 

Combining this with Lemma |5. 11 Fiirstenberg's Theorem gives the required assertion. □ 

Proof of Theorem \3.2\ Again by Fiirstenberg's Theorem and Lemma 15.11 proving that 
7(A) > for all A € C for the case that the support of contains two points of T that 
are not diametrically opposed, is reduced to checking condition (|20p. Each element of the 

projective space P(C^) is of the form v with v = ox v = 

In terms of x, z, introduced above, the condition that |a — 6| ^ {0,7r} can be written as 
xz^{-l,l}. 



, for some a £ 
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Case I: Let v = ■ '^^^ action of the operator E, defined in p!7|). on v has 

/L(xz-1)\ _ f-((xz)'^-l)\ 

the direction of t , while E'^v has the direction of t . Thus 

VI/ V 1 y 

{I,E,E^} C Gx^^ maps v into three different elements in P(C^). 

Case II: Let v = ' '^ith a E C. Acting on v with the operator D from (jl6p results 

f ^ \ ( ^- 

in the direction of r _ _ | , while has the direction of | o^xz 

I -1 1 — xz) + axz 



,t / \l + ar/t(l — xz) 

T _ _ 

Defining the map F such that F : c ^ -{\ — xz) + cxz, one sees that F has a single 

fixed point at c = r/t. Since the second iteration : c i— > -(1 — (xz)^) + c^xz)"^, has the 

same value r/t as its only fixed point, we deduce that {c, F{c), F^{c)} are pairwise different 
except when c = r/t. Thus iterations of the operator D take v into at least three different 
directions, unless a = r/t. 

cxz 

On the other hand, the map H : c ^ r- — - has fixed points 0, — t/r, which 

1 + cr/t{\ — xz) 

are also the fixed points of : c i— > , In particular, {/, EjE"^} C G\ „ 

1 + cr /t[l — [xz)"^) 

map the direction vector ^ to three different elements in P(C^). This proves the 

required condition for strong irreducibility of G^.^ and once more the result of Lemma |5. II 
and Fiirstenberg's Theorem finish the proof. □ 

6 A Unitary Dimer Model 

In this section we study a unitary version of the Dimer model, which is obtained from the 
Anderson model ([6]) by doubling up the random phases. More rigorously, let a probability 
measure ^ on T be given, define (7 : T ^ as g{6) = {0, 9) and let p. be the probability 
measure supported on the diagonal of induced by /i through g: fl{B) = fi{g~^{B)) for 
Borel sets B inT'^. We introduce the probability space {Cl, T, IP), where is identified with 
(T^)^, T is the cr-algebra generated by cylinders of Borel sets in and IP = (S)fcez A- 

For Lo = {uik)kez £ ^, the random phases 9'^, n G Z, used in ^ and ([6]) to define Uuj 
are now chosen as 

(^2/0 (^2k+i) = ^k, k e z. 

For this model we again prove that the almost sure absolutely continuous spectrum is 
empty. We also show that the case of a Bernoulli measure 

fJ,=pSa+q5b, p + q = l, a,beT, a ^ b, (21) 

gives rise to additional critical quasi-energies, as this is the case of least randomness. 

The following theorem states that for any non-trivial distribution on T, i.e. {a, 6} C 
supp fj, foi a b, the Lyapunov exponent is positive for all but a finite set of quasi-energies, 
given by 

M := {-a,-b}u{ManMb), 
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where 

Ma = {arccos(r^) — a, 27r — arccos(r^) — a, arccos(r^ — t^) — a, 27r — arccos(r^ — t^) — a}, 
Mb = {arccos(r^) — 27r — arccos(r^) — 6, arccos(r^ — t^) — 6, 27r — arccos(r^ — t^) — h}. 

An immediate consequence is that the almost sure absolutely continuous spectrum of these 
operators is trivial. 

Theorem 6.1. If {a,b} C suppfi, then for all A G T\M, the Lyapunov exponent 7(A) is 
strictly positive. In particular, S^c = 0- 

As before, we will use Fiirstenberg's Theorem [2] to prove positivity of Lyapunov expo- 
nents. Let Ga,^ be the closed group corresponding to Gx,^ from the previous section. 

We will show that Gx^^ is both non-compact and strongly irreducible for all A outside 
of M. As both of these properties carry over to larger groups, we may assume for the 
rest of the proof of Theorem 16.11 that supp /i = {a, b}. Thus Gx,^ is generated by just two 
matrices, T{9, 9) and T(r/, rj), where 9 = a + X, rj = b + X. 

In order to prove Theorem 16. H we start by mapping the problem into a somewhat 
simpler form. In order to simplify the notation, we again let x := e~*^, z := e~*'', and let 

2r^ 1 

p, pi be the two eigenvalues of T(9, 9). Since tr T{9, 9) = ^(x -l-x), and det Ti9, 9) = 

t^ 

detr(?7, 77) = 1, we have the following cases 

p = p]-, |p| = l, when |trr(6',6l)| < 2 

p = — > 1, when |trr(6',6l)| > 2 (22) 
Pi 

p = pi,p2^1, when |trr(6',6')| = 2. 



This allows us to introduce the transformation N given by 
N 




x + p 
r 

^2 

Using that ix+ p){x+ pi) = — ^rfl — x)^, we deduce that det = {x+p){pi— p). Therefore, 
is invertible as long as |trr(^,^)| 7^ 2. Moreover, 

E = NT{9,9)N-^=(P^ (23) 

A short calculation shows that the elements of = NT{ri, ri)N^^ are given by 

^11 = t^{p-l/p) [^"^^^^ + " + ^""^ " ''^^ + 

'2,77' 

F12 = F2I = ^3(p_i/^) [^(^) - ^(^) + ^(^^)]' 

^22 = ,2. x [2r-'(l + -) - {zx + zx) - -{z + z)]. (24) 

t^p-l/p) p p 

Notice that since 9 ^ tj, F12 = if and only if either rj = or 9 = 0. 

Since proving non-compactness and strong irreducibility of Gx,^ is equivalent to proving 
the same properties for the group Hx,^ generated by the matrices E, F, we will use the 
latter, somewhat simpler matrices whenever it helps simplifying the proofs. 
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Lemma 6.1. For all A € T\{— a, —b}, the group G\^^ is non-compact. 

Proof. Since A 7^ —a, we have that ^ 7^ and thus trT{6, 9) 7^ —2. Therefore, the preceding 
discussion suggests the proof should be divided into the following cases; 

Case I: tr T{e, 6) = 2. By ([221) we have p = pi = 1 and since by definition T{e, 0) I 
it follows that there exists a non-singular matrix R such that 



RT{9,e)R 



-1 



P 1 
p 



Since \\[RT{6,6)R^^]'^\\ grows with n, the group generated by RT{6,6)R^^ , RT{r],r])R^^ 
is non-compact, which implies that Gx^^j^ is non-compact. 

Case II: \tr T{6, 6)\ > 2, again by (I22p . T{6, 6) has an eigenvalue p > 1 which gives the 
required result. 

Case III: |trr(6',6')| < 2. In this case equations ([22]), ([23]) give 



E 



e'y 
e-'y 



,r2 1 



with y = arccosf-^ 7t(x + x)), and y G (0, vr). Equations lead to 

2t^ 

^ a>0,/?EM,cGT. 

Using that detF = 1, it follows that a > 0. 

Now we follow a strategy outlined in [10] to show that there exists a sequence of elements 



in Hx^^ with unbounded norms. In order to do so, we note that any element of P(C^) can 



be written in the form 



e(n,^,) = y^-,u j , (^^, G [0, vr) X [0, ^). 

Therefore, for any element e{u^v) of PCC'^) we have 



ae 



cos(t;) + /3e-^" sin(t>) 



V«e"'^''^"^ sin(z;) + /3e^" cos(7;) y ' 
Using that — = 1, we get 

1 Ip — 1 = 2/?^ + 4a/3 cos(2ti + c) cos(t;) sin(u). 



^ r Sjfx) — ^(z) + 9(zx) ..... 1-1 ,1 , , 

Smce p = -5- -— ^ with distmct x,z and neither equals 1 under the 

Q{p) 

current assumptions, we see that /? 7^ 0. 

■V il I I 

f2 1 ^ ^2 : 



In the case /? > 0: If cos(?;) sin(t!) = 0, then ||-Fe(„ |p - 1 > for all u G [0, vr). While 
for cos{v) sin(ti) > 0, ||-Fe(u„)||^ — 1 > is equivalent to 



cos(2u + c) > 



4a cos(^;) sin(u) 
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—13 

In particular, the condition cos(2n+c) > - — guarantees that I li^ef^ ^^-i I P > Defining 

4a ^ ' ' 

Kj^ := {u G [0,7r) : cos(2n + c) > -r^}, we see that ji^+| > vr/2 and for ah u € we 

4a 

have ||-Fe(„^„)|p > 1 + 

Similarly, for cos(t>) sin(t)) < 0, let := {u £ [0, tt) : cos(2ti + c) < -|— }. Then for all 

4a 

u G we have ||Fe(„_^)|p > 1 + and |i^_| > 7r/2. 



Hence, given any v € [0, vr), there exists an interval Ky C [0, vr), i.e. an interval in ^ 
such that \Ky \ > 7r/2 and ||Fe(„ t,)||^ > 1 + /3^, for all u G K^. Therefore, starting with an 
appropriately chosen vector e(„ ,;) G P(C^) such that | |-Fe(u | p > 1 + applying F will 
result in a vector ce(^ui,vi) with c > 1. Now we apply E as many times as required to get a 
vector ce(s_^^) (or — ce(s,;^)) such that u G i^t,^. Iterating this process gives a sequence of 
vectors with unbounded norms. The case /? < is treated similarly. Thus we have proved 
the non-compactness of Hx,i_i, and consequently that of G^,^. □ 

The next step is proving that Ga,/^ is strongly irreducible for all A outside the set M. 

Lemma 6.2. G^.^ is strongly irreducible, for all A G T\M. 

Proof. Since we already proved that G^,^ is non-compact (Lemma 16. ip . it suffices to show 
that for all v G P{C^),i^{gv : g G Gx,^,} > 3. 

We first note that = 1 implies that \tiT(6,9)\ G {0,2}, which in turn implies that 
A G {—a, —b} U Ma- Therefore, the condition A ^ {—a, —6} U Ma gives that ^ 1. Hence, 
{I,E,E^} C Ga,^ maps every v G P(C^) to three different directions unless v coincides 

with either ( ^ ) or (^] ■ At this point, we note that the diagonal elements of the matrix F, 



given in (j24p . vanish simultaneously only if A G {—a, —6}. Thus, iterations of the operator 
F, followed if necessary with iterations of E, maps each of the latter directions to at least 
three different elements of P(C'^). This proves strong irreducibility of the group Hx ^j^, for 
all A ^ {—a, —b} U Ma, which gives the corresponding result for Ga,/^. A similar argument, 
replacing the rules of a and b, gives a similar assertion for A ^ {—a, — 6}UMb, thus finishing 
the proof. □ 

Proof of Theorem \6.1[ The previous two lemmas combined with Fiirstenberg's Theo- 
rem, immediately give 7(A) > for all A G T\M. Using Theorem 12. H we deduce that 

Sac = 0. □ 

Even though the finiteness of the set M is more than enough to prove the absence of 
absolutely continuous spectrum for the unitary dimer model, M is by no means optimal. 
Determining whether or not a certain element of M is, in fact, a critical quasi-energy of U^j 
requires further analysis. Nevertheless, the proof suggests that if the support of /i contains 
three or more points than generically M = 0. Even for the Bernoulli unitary dimer model 
where supp/i = {a, b} the analysis is likely to fall into a number of different sub-cases. 
However, for generic choices of 0,6 we see that Ma Mj, = 0, thus M = {—a,—b} and 
indeed the situation where A G {—a, —b} is readily accessible. Guided by the proof of 
Theorem 2.2(i) of [11] for the self-adjoint dimer model we prove that, 

Proposition 6.1. For a probability measure /x given by h21\). we have the following 
(i) // |a - 6| G a{S), then 7(-a) = 7(-5) = 0. 
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(ii) // \a-b\ e p{S), then both 7(-a) > and 7(-6) > 0. 

Proof. For A = —a, we have that T{9,9) = —I, and r] = b — a. Since tTT{r],rj) = 
'2t'^ 2 cos(/y) 

p , then |trT(r/, 'q)\ <2 when b — a € cr{S) and |trT(r/, r])\ > 2 if 6 — a G p(«S'). 

Now, let rUn ■= : 1 < A; < n, ^gfc = then P-almost surely 

lim = q. 

n^oo n 

This along with the fact that 

lim 1 1 [T(r/, 77)]"^" 1 1 1/™" = max jr^l, 

where rj are the eigenvalues of T(r],rj), gives the results of the proposition. □ 

7 Continuity of the Lyapunov Exponent 

In this section we prove that, away from the critical points, the Lyapunov exponent is a 
continuous function of the spectral parameter A. The proof of this fact is similar to the 
one given in ^ for the self-adjoint case. 

First, for a compact interval I of quasi-energies with positive Lyapunov exponents, we 
define the function 

$(A,t;) =E(ln^^), E P(C2), A G /. 

\\v\\ 

Where Ta := r(0^„(A), 

^2n+iW) denotes the transfer matrix defined in ([7]) with the depen- 
dance on uj being suppressed in order to simplify the notation. The next lemma establishes 
a couple of properties of $(A,-u). 

Lemma 7.1. (i) The mapping v 1— > $(A,v) is continuous on P(C^). 
(ii) There exists a constant C such that, 

sup \^{\,v)-^{Xi,v)\<C\e'^-e'^^, A, Ai G /. 

i)GP(C2) 

Proof, (i) From ([7]) one sees that the norm of Tx is uniformly bounded for all X,uj. Con- 
sequently, we also have a uniform bound on ^(A,^;). The assertion is then obtained using 
the dominated convergence theorem. 

(ii) First we note that for all v G P(C^), 

|cl>(A,^)-$(Ai,^)|<E(|lni^|). 

Since IdetTxT^^^I = 1, it follows that 

MX,v) - ^Xi,v)\ <E{ln\\TxT-^'\\). (25) 
On the other hand, one has 

||TAr,-ii|<||r,-/||||r,-r,-i|| + i. 
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Since all norms on GL(2,C) are equivalent, there exists a constant Ci such that 

||TAr,-i||<Ci||rA-r,-;i||^ + i, (26) 

where ||j4|||, = Ylf j=i\^ij\'^ denotes the Frobenius norm of the matrix A = {aijjfj^^. 
Prom dZD, it is easy to see that \\Tx - T~^^\\f < C2\e'^ - e*^i|. Combining this with 1^ 
and (f25l) gives the required result. □ 

Before proving the main result of this section, we recall a general fact: If fix denotes 
the probability measure on GL{2,C) induced by Tx and Gx,^ is non-compact and strongly 
irreducible for all A G /, then there exists a unique distribution vx on P(C^) that is invariant 
with respect to fix- A proof of this fact can be found in [2]. Moreover, we have 

Lemma 7.2. For A € /, the mapping A i— > is weakly continuous. 

Proof. We start by showing that if the sequence {Xn}nez converges to A, then the corre- 
sponding measures fiXn converge weakly to fix- First recall that for all cj € fi, ||Tx„(u;) — 
Tx{io)\\ < dn, where d„ = C|e*^» - e*^| for some C G M. Now, let B C GL(2, C) such that 
the boundary of B has zero measure with respect to fix, i.e. fix{dB) = 0. For such a set 
we have that 

Ifix^B) - fixm = |P[r-i(i?)] - nTx\B)]\ 

< F[T-^HB) n {T-\B)y]+F[T-\B) n {T^^}iB)y] 

It is not difficult to see that 

nTx„\B) n {T^\B)y] < F[{u; G Q : Tx{io) G B'',dist{Tx{uj),dB) < 

< f-ix[{A : dist{A,dB) < dn}] 

Taking the limit as n — > oo and using dominated convergence one sees that 

hm P[T-i(i?) n (,T-\B)r] = fixidB) = 0. 

n—too " 

Using a similar argument one gets that ¥[T^^{B) n {T^^ {B)y] — > as n — > oo. Therefore, 
we have 

hm fix„{B) = fix{B). 

Since this is true for any set B with fix{dB) = 0, weak convergence of to follows pO] . 
In order to get the weak convergence of z^a„; we use the fact that the set of invariant 
measures on P(C^) is compact in the weak* topology [20], thus every subsequence of {i^A„} 
has a weakly convergent subsequence and since the limit of each of those subsequences is 
invariant with respect to fix it equals ux by uniqueness of the latter. A short contradiction 
argument shows that vxn has to converge weakly to ux- □ 

The Lyapunov exponent can be expressed in terms of the mapping $ and the measure 
fx as 

7(A) = j ^{\,v)dux{v). (27) 

Now we are ready to prove that for any interval / for which Ga,^ is non-compact and 
strongly irreducible for all A G /, we have 
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Theorem 7.1. The Lyapunov exponent 7(A) is a continuous function of X £ I. 



Proof. Let As/ and choose {An} C I a sequence of quasi-energies such that A; 
n — > CX3. It follows that 



A as 



lim 7(A„) = lim / ^{Xn,v)diyx„{v) 



00 n— >oo 




= lim [/ ^X,v)dux„iv)+ mXn,v)-^X,v))di^xAv)] 

n^oo J J 

= / ^{X,v)dux{v), 



with the last equality following from Lemma 17.21 along with part (ii) of Lemma I7.1[ □ 
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